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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 
TECHNICAL  NOTE  3016 


ANALYSIS  OF  TURBULENT  HEAT  TRANSFER  AND  FLOW  IN  THE 
ENTRANCE  REGIONS  OF  SMOOTH  PASSAGES 
By  Robert  G.  Deissler 


SUMMARY 

A previous  analysis  for  fully  developed  turbulent  heat  transfer  and 
flow  with  variable  fluid  properties  is  extended  and  applied  to  the 
entrance  regions  of  smooth  tubes  and  parallel  flat  plates.  Integral 
heat-transfer  and  momentum  equations  are  used  for  calculating  the  thick- 
nesses of  the  thermal  and  flow  boundary  layers.  The  effect  of  variable 
properties  is  determined  for  the  case  of  uniform  heat  flux,  uniform 
initial  temperature  distribution,  and  fully  developed  velocity  distri- 
bution. A number  of  other  cases  in  which  the  fluid  properties  are  con- 
sidered constant  are  analyzed.  The  predicted  Nusselt  numbers  for  air 
with  a uniform  wall  temperature  and  uniform  initial  temperature  and 
velocity  distributions  agree  closely  with  experimentally  determined 
values. 


INTRODUCTION 

Flow  and  heat  transfer  in  the  entrance  regions  of  channels  have 
been  subjects  of  considerable  Interest.  The  laminar  flow  case  has 
been  analyzed  in  references  1 (pp.  299-310),  2,  and  3 (pp.  451-464), 
for  instance,  and  the  turbulent  case,  in  references  4 and  5.  In  these 
analyses  constant  fluid  properties  war's  assumed,  and  uniform  wall 
temperature  was  postulated  in  the  cases  in  which  heat  transfer  was  con- 
sidered. In  the  turbulent  flow  cases  the  flow  was  assumed  to  be  tur- 
bulent at  all  points  along  the  passage. 

The  analysis  in  reference  4 is  for  a Prandtl  number  of  1 and  is 
based  on  an  assumed  l/7-power  velocity  profile  and  the  Blasius  resis- 
tance formula.  The  relative  changes  in  heat-transfer  coefficient  along 
the  tube  predicted  in  that  analysis  are  approximately  correct,  but  the 
absolute  values  of  the  heat-transfer  coefficient  at  the  end  of  the 
entrance  region  are  about  20  percent  too  low.  The  analysis  given  in 
reference  5 is  applicable  to  cases  in  which  the  heat  transferred  by 
turbulence  can  be  neglected,  as  might  occur  at  very  low  Prandtl  numbers. 
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The  analysis  given  herein,  which  was  made  at  the  NACA  LewlB  labora- 
tory, is  an  extension  of  the  analysis  given  in  reference  6 for  fully 
developed  heat  transfer  and  flow  with  variable  fluid  properties  to 
the  entrance  regions  of  smooth  passages.  In  calculating  the  growth  of 
the  thermal  and  flow  boundary  layers  in  the  entrance  region,  integral 
equations  for  heat  transfer  and  momentum  are  used . For  obtaining  the 
velocity  and  temperature  distributions  in  the  boundary  layers,  the  same 
assumptions  are  made  for  solving  the  turbulent  transfer  equations 
as  were  made  in  reference  6 . 

The  flow  is  assumed  to  be  turbulent  at  all  points  along  the 
passage,  as  in  previous  analyses,  so  that  the  analysis  should  be 
applicable  when  disturbances  occurring  at  the  entrance  are  sufficient 
to  produce  a turbulent  boundary  layer.  According  to  Prandtl's  assump- 
tion, which  is  in  reasonable  agreement  with  experiment,  in  the  presence 
of  a laminar  boundary  layer  near  the  entrance,  the  turbulent  portion  of 
the  boundary  layer  behaves  as  th' ugh  the  boundary  layer  were  turbulent 
all  the  way  from  the  entrance  (ref.  7,  p.  74).  The  present  calculations 
may  therefore  be  applicable  to  the  turbulent  portion  of  the  boundary 
layer  even  when  a laminar  boundary  layor  exists  near  the  entrance. 

The  cases  considered  in  the  present  analysis  are  given  in  the 
following  table: 


Phenomenon 

Wall 

boundary 

condition 

Initial 

temperature 

distribution 

Initial 

velocity 

distribution 

Passage 

Prandtl 

number 

Proper- 
t ies 

Heat 

transfer 

Uniform 
heat  flux 

Uniform 

Fully 

developed 

Tube 

0.73 

Variable 

Heat 

transfer 

Uniform 
heat  flux 

Uniform 

Fully 

developed 

Parallel 

flat 

plates 

.73 

Variable 

Heat 

transfer 

Uniform 
wall  tem- 
perature 

Uniform 

Fully 

developed 

Tube 

.73 

Constant 

Heat 

transfer 

Uniform 
heat  flux 

Uniform 

Uniform 

Tube 

.73 

Constant 

Heat 

transfer 

Uniform 
heat  flux 

Uniform 

Uniform 

Parallel 

flat 

plates 

.73 

Constant 

Heat 

transfer 

Uniform 
wall  tem- 
perature 

Uniform 

Uniform 

Tube 

.73 

Constant 

Heat 

transfer 

Friction 

Friction 

Uniform 
heat  flux 

Uniform 

Fully 

developed 

Uniform 

Uniform 

Tube 

Tube 

Parallel 

flat 

plates 

.01 

Constant 

Constant 

Constant 

CZ-1  back: 
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Although  moat  of  the  cases  In  the  preceding  table  are  for  constant  fluid 
properties,  the  greator  portion  of  the  report  Is  concerned  with  the 
first  two  cases,  in  which  the  properties  are  variable. 


ANALYSIS 

For  calculating  host,  transfer  and  friction  In  the  entrance  region 
go  of  ducts,  the  usual  boundary  layer  assumptions  are  made;  that  Is,  It  Is 

cn  assumed  that  the  effects  of  heat  transfer  and  friction  are  confined  to 

fluid  layers  close  to  the  surface  (thermal  and  flow  boundary  layers, 
respectively).  The  temperature  and  velocity  distributions  outside  the 
boundary  layers  are  assumed  uniform,  and  the  total  temperature  and  total 
pressure  are  constant  along  the  length  of  tho  duct  for  the  region  outside 
the  boundary  layers.  More  exact  analyses  (ref.  2)  indicate  that  these 
assumptions  are  valid,  even  for  laminar  flow,  except  In  the  region  at  a 
distance  from  tho  entrance  where  the  boundary  layer  fills  a large  portion 
of  the  tube.  For  that  region,  however , the  Nusselt  numbers  and  friction 
factors_have  values  very  close  to  the  values  for  fully  developed  flow. 

In  tho  following  analysis,  flow  in  round  tubes  and  between  parallel 
flat  plates  will  be  considered.  The  effect  of  frictional  heating  will 
be  neglected  and  then  later  investigated  to  3ome  extent  in  the  section 
Effect  of  frictional  heating  on  entrance  length. 


* 


Velocity  and  Temperature  Distributions  in  Boundary  Layers  for  Gases 


For  obtaining  the  velocities  and  temperatures  in  the  flow  and 
thermal  boundary  layers,  tho  differential  equations  for  shear  stress 
and  heat  transfer  can  be  written  as  follows: 


„ du  „ du 

— p.  j j 

dy  dy 


- k | - PCPe£*  | 


(The  symbols  U3ed  in  this  report  are  defined  in  appendix  A.)  In  the 
preceding  equations  t and  are  the  eddy  diffusivities  for  momentum 

and  heat  transfer,  respectively,  the  values  for  which  are  dependent  on 
the  amount  and  kind— of  turbulent  mixing  at  a point.  When  written  in 
dimensionless  form,  these  equations  become 

x p du+  P e du+ 

T0  Po  dy+  + p0  Po/pO  dy+ 


r 
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and 

4_  = fjL  _L  P_  fj2_  ch  \ dt+ 

*0  ~ \kO  + p0  CP0  ^o/poJ  dy+ 

where  the  subscript  0 refers  to  values  at  the  wall. 

Assumptions . - The  following  assumptions  are  made  in  the  use  of  the 
previous  equations  for  obtaining  velocity  and  temperature  distributions 
in  the  flow  and  thermal  boundary  layers: 

1.  The  eddy  diffusivities  for  momentum  e and  heat  transfer 

are  equal.  Previous  analyses  for  fully  developed  flow  in  tubes  based  on 
this  assumption  yielded  heat-transfer  coefficients  and  friction  factors 
that  agree  with  experiment  for  Reynolds  numbers  above  15,000  (ref.  6). 

At  low  Reynolds  or  Peclet  .numbers  (Pe  = RexPr),  the  ratio  of  eddy 
diffusivities  appears  to  be  a function  of  Peclet  number-. (ref  8) ; but 
for  Reynolds  numbers  above  15,000,  it  is  nearly  constant  for  gaseB,  at 
least  for  fully  developed  boundary  layers.  The  same  assumption  is  made 
here  for  developing  boundary  layers. 

2.  The  eddy  diffusivity  t is  given  by 

e = n2uy 

in  the  region  close  to  the  wall  (y+<26),  and  by  the  Karman  relation 

_ 2 (du/dy)3 * * 6 

C A 2 

(d2u/dy2) 

in  the  region  at  a distance  from  the  wall  (y+>26).  These  expressions 
have  been  experimentally' verified  for  fully  developed  boundary  layers 
with  variable  properties  in  reference  6.  It  is  assumed  here  that  they 
apply  also  to  developing  boundary  layers.  The  quantities  n and  x are 
experimental  constants  “having  the  values  0.109  and  0.36,  respectively. 

3.  The  variations  across  the  flow  and  thermal  boundary  layers  of 
the  shear  stress  t and  the  heat  transfer  per  unit  area  q have  a 
negligible  effect  on  the  velocity  and  temperature  distributions.  It  is 

shown  in  figure  12  of  reference  6 that  the  assumption  of  a linear  varia- 
tion of  shear  stress  and  heat  transfer  across  the  boundary  layers 

(t  or  q = 0 at  the  edge  of  the  thermal  or  flow  boundary  layer)  gives 

very  nearly  the  same  velocity  and  temperature  profiles  as  those  obtained 
by  assuming  uniform  shear  stress  and  heat  transfer  across  the  boundary 
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layers  for  values  of  b+  or  between  500  and  5000. For  small 

values  of  8^*,  such  as  occur  very  near  the  entrance,  the  effect  of 

variable  heat  transfer  will  be  checked  in  appendix  B and  In  the  section 
Effect  of  variation  of  heat  transfer  across  boundary  layer  on  Nusselt 
number  for  small  5j1+  and  constant  fluid  properties. 

4.  The  molecular  shear  stress  and  heat-transfer  terms  in  the 
equations  can  be  neglected  In  the  region  at  a distance  from  the  wall 
(ref.  6,  fig.  14). 

5.  The  static  pressure  can  be  considered  constant  across  the 
passage . 

6.  The  Prandtl  number  and  specific  heat  can  be  considered  constant 
with  temperature  variation.  The  variations  with  temperature  of  the 
specific  heat  and  Prandtl  number  of  gases  are  of  a lower  order  of  magni- 
tude than  the  variations  of  viscosity,  thermal  conductivity,  and 
density.  The  variations  of  viscosity  and  thermal  conductivity  with 
temperature  are  assumed  to  be  given  by  the  equations 

p/pQ  * k/k0  = (t/t0)0-68 

From  the  definitions  of  the  quantities  involved  it  can  be  shown  that 

t/tQ  - 1 - pt+ 

Solution  of  equations.  - The  details  of  the  solution  of  the  equa- 
tions for  shear  stress  and  heat  transfer  under  the  foregoing  assumptions 
are  given  in  reference  6.  The  resulting  relations  between  t+  and  y+ 
and  u+  and  y+  are  given  in  figures  1 and  2.  Positive  values  of  P 
correspond  to  heat  addition  to  the  gas;  negative  values,  to  heat 
extraction.  For  the  cases  in  which  the  fluid  properties  are  considered 
constant,  the  curves  for  P = 0 are  used.1"  The  plotted  values  of  t+ 
and  u+  are  to  be  used  in  the  thermal  and  flow  boundary  layers,  that 
is,  for  values  of  y+<  &h+  or  y+<  & + - For  values  of  y+  greater 


^Figure  12  of  reference  6 applies  to  developing  boundary  layers  if 
the  symbols  Tq+  in  the  figure  are  replaced  by  8+  or  8^+. 

^The  case  in  which  P = 0 is  a limiting  one  which  can  be  approached 
as  closely  as  desired  by  making  the  temperature  difference  small.  In 
the  limit,  where  the  temperature  difference  is  zero  or  infinitesimal,  the 
properties  must  be  constant. 


6 


NACA  IN  3016 


than  these  values,  t+  and  u+  are  constant  and  have  the  values  at  &^+ 
and  b+ . The  relations  between  5^+  and  X/D  and  &+  and  X/D  will 
be  obtained  in  the  following  sections. 


Development  of  Thermal  Boundary  Layers 


Heat-flow  equations  for  thermal  boundary  layers.  - The  heat -flow 
equation  for  the  thermal  boundary  layer  in  a round  tube  can  be  obtained 
from  the  following  diagram: 


Energy  flows  into  the  differential  annulus  by  convection  through  plane  1 
and  leaves  through  plane  2.  In  addition  to  this  energy  there  is  a radial 
flow  of  energy  (heat)  from  the  tube  surface  by  conduction  and  a radial 
flow  by  convection  at  b^.  No  radial  heat  flow  due  to  temperature 

(gradient  -takes  place  at  because  the  temperature  gradient  is,  by 

definition  of  the  thermal  boundary  layer,  zero  at  the  edge  of  the 
boundary  layer.' 


Equating  the  heat  energy  entering  the  annulus  to  that  leaving 
gives,  for  constant  . Cp, 


q0  2rt  Tq  dX  + 


pu  2nr  dy 


tPu  Cp  g 2itr  dy ) + 


(f 


pu  2rtr  dy 


cp  8 t5 


’h 


tpu  Cp  g 2«r  dy 


> 


« 
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where  the  subscripts  on  the  parentheses  refer  to  planes  1 and  2 in  the 

diagram,  or 


4.,  r(,  dX  = d 


-'p  ^'6 


cp  c tpu(rQ  - y )dy 


Pu(r  - y)d 


where  the  different  in is  of  the  integral--  indicate  changes  in  the  X-direction. 
Equation  (])  iw.  the  Jo.- i red  equation  for  flow  in  tuber  relating  the  thermal- 
t> 'lundary- layer  thiekne.  s to  the  distance  along  the  tube. 

Ait  equation  corresponding  to  equation  (i)  can  be  easily  derived  for 
flow  between  fiat  plates.  Tr.is  equation  is 


0,-  dX  - d 


(/ 


cp  g tpu  d.vJ  - cp  tb  g d 


Pu  dy 


Equation  (2)  in  the  incompressible  form  has  been  used  extensively 
for  calculating  the  growth  of  the  thermal  boundary  layer  on  a flat 
plate  (r-f.  7,  p.  Sfj. 

Uniform  heat  flux,  variable  fluid  properties  (gases).  - The  total 
temperature  outside  the  thermal  boundary  layer  does  not  vary  with  X 
inasmuch  as  no  heat  penetrates  the  region  outside  the  boundary  layer. 
(The  differences  between  total  and  static  temperature  are  small  except 
at  high  Mach  numbers.)  For  uniform  heat  flux  at  the  wall,  equation  (l) 
can  then  bo  integrated  to  give 


*0  rG  x 


(t  - tb)  Pu(r0  - y)dy 


Equation  (3)  can  be  written  in  dimensionless  form  as 


(tb+  - t+)  ^ u+(*o+  ’ y+)dy+ 


r 


'a 
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where,  from  the  definitions  of  the  variouo  quantities  involved  and  the 
perfect  gas  law  (constant  pressure ) , 

p/p0  --  1/(1  - pt+)  (5) 


When  equation  (4)  is  to  be  used,  values  of  r,j+  and  3 must  first  he 

fixed.  Values  of  X/D  can  then  he  calculated  for  various  values  of 
5j1+.  The  relations  between  t 4 and  y+  and  u+  and  y+  are  obtained 

from  figures  1 and  2. 

The  equation  for  fiat  plates  corresponding  to  equation  (4)  is 
olt Mined  from  equation  (2)  and  is  found  to  be 

(tb+  - t+)  u+  dy+  (6) 


A 

D 


tr 


where  D is  the  equivalent,  diameter  f.r  flow  between  flat  plates,  or 
twice  the  plate  spacing. 

Uniform  wall  temperature,  constant  fluid  properties.  - For  flow  in 
tubes,  equation  (l)  can  bo  written  as 


where  the  quantity  in  brackets  is  also  the  upper  limit  on  the  integral. 
For  constant  fluid  properties  and  uniform  wall  temperature,  equation  (7) 
can  be  written  in  dir.onsi  nless  fora  as 


1 — 1 
1 1 

T P bv. 

Jo  ro+p  d 

— 7 / (tb+  - t + )p  u+(r0+  - y+)  dy+ 

ro  Jo 

(8) 


where  the  shear  stress  xf]  is  allowed  to  vary  with  X.  But 

tb/to  = ti/tu  - 1 ' P 


since 


t{j  does  not  vary  with  X.  Therefore, 

0 = d/VJd-^/tQ) 
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.1 


ro 

co 

■%J 


Substituting  this  expression  for  p into  equation  (8)  results  in 


Equation  (y)  gives  the  relation  letwoon  i./V>  and  b^+  for  round  tubes 

with  uniform  wall  temperature  and  constant  fluid  properties.  Its  use 
is  similar  to  that  of  equation  (4). 


Nuaselt  numbers  and  Reynolds  numbers.  - It  can  easily  be  shown 
(see  ref.  i>)  from  the  definitions  of  the  quantities  involved  that  the 
Nussult  number  and  Reynolds  number  (properties  variable  or  constant) 
for  circular  tubes  with  the  fluid  properties  evaluated  at  the  wall  tem- 
perature are  given  by 

Pr 


Nu. 


(10) 


■and 


Jub  r 


(11) 


whore 


\ y»  . * , + / -f  + \ 

r0  t u (rL,  -y  ) 
1 - pt + 


dy+ 


'r,-,  +/  + +\ 

u (r0  ~>v  ' 

i-pt  + 


d.v+ 


(12) 


and 


r0 


u+(r0+-y+)dy+ 
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wble'.\  by  combination  with  equation  (12),  becomes 


Ub+  - ~77  (1-Ptb  + ) 


nb 


h + / + + \ 

u (ro  -y  ) , + i 

dy  + 


1-pt 


1-Pt& 


Dh 


u+(rn+-y+)dy+ 


(13) 


The  foregoing  equations  uru  for  a circular  tube.  The  corresponding 
equations  for  flow  between  flat  plates  are 


When  values  of  Nu^  and  Boq  are  calculated,  values  of  rQ+  and  P 
must  be  fixed  as  in  equation  (-1).  Values  of  Nu0  find  Rgq  are  thus 
obtained  for  given  values  of  X/D  and  p.u 


Unfortunately,  p varies  along  the  passage.  It  will  later  be  seen 
that  qj_ ' is  more  convenient  for  use  as  a heat -flux  parameter. 


CZ-2  back 


NACA  TN  3016 


11 


ro 

OP 

C.n 

-j 


The  temperature  and  velocity  distributions  from  figures  1 and  2 can 
be  used  for  these  quantities  except  for  the  region  outside  the  thermal 
boundary  layer.  Thu  temperature  distribution  is,  of  course,  uniform 
outside  the  thermal  boundary  layer  (t+  = to+).  When  the  velocity  distri- 
bution in  that  region  is  calculated  for  the  case  in  which  the  velocity 
distribution  is  fully  developed,  it  should  be  remembered  that  the 
properties  are  constant  and  have  the  values  at  &j1+.  For  values  of 
S>b+>26,  the  velocity  distribution  curves  outside  the  thermal  boundary 

on  the  semi logarithmic  plot  are  straight  lines  having  the  slope  of  the 
curves  at  6b+. 

Nunselt  numbers  and  Reynolds  numbers  with  properties  based  on  tem- 
peratures other  than  the  wail  temperature  can  be  obtained  by  use  of  the 
relations 


tb/t0  = 1-Ptb+ 

H/to  = 

Pb/P0  - i/d-Ptb+) 

MbA,;  - kpA'o  = (l-iitb  + )d 

and  so  forth,  whore  the  perfect  gas  law  with  constant  static  pressure 
across  the  tube,  constant  Prandtl  number,  and  constant  specific  heat 
have  been  assumed  as  in  reference  6.  The  exponent  d is  obtained  from 
viscosity  data  and  is  found  to  be  about  0.68  for  air  and  most  gases. 

The  following  relations  for  Reynolds  numbers  then  result: 


Reb 


Re 


= 


1 * Rob 


1 

(tbA0)d+1 

(tb/to)d 

(ti/to)d 


(18) 

(19) 


where  the  continuity  relation  Pb  ub  = Pb  ub  has  been  used.  Similarly, 


Nu, 


Nub  -- 


Ab/^o) 


(20) 
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and 


Nujl 


^lAb^ 


(21) 


Another  useful  parameter  is  ' , which  can  be  calculated  from 

P(tb/t0) 

qi  = — : 

% (ti/tQ) 


(22) 


Equation  (22)  can  bu  obtained  from  the  definition  of  q1 1 and  the 
perfect  gas  law  for  constant  pressure,  t0/tk  = Pb/pQ>  for  P-^  u^  is 
constant  along  the  tube.  The  important  property  of  q^  1 is  that 
it  does  not  vary  along  the  length  of  the  tube. 


It  is  desired  to  calculate  the  variation  of  bulk  temperature  and  wall 
temperature  along  the  tube  for  uniform  heat  flucc.  It  can  be  shown,  by  use 
of  the  expressions  q0  - h(tL)-tb)  and  nDXq0  = Pg(n/4)D2  ut  cp(tt)-t1), 

that  dimensionless  bulk  temperature  and  wall  temperature  parameters  are 
given  by 


(23) 


(24) 


Equations  (23)  and  (24)  can  be  used  for  flow  between  flat  plates  as 
well  as  for  tubes. 

The  average  Nusselt  number  for  a given  X/D  i3  calculated  from 


Kui,av  ~ px/D  (2^) 

/ a (x/d) 

Jo  Nui 


This  equation  is  consistent  with  the  definition 
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(to  - tb)av 


(t0-tb)dx 


X 


(26) 


which  ia 
and  bulk 
equation 


the  usual  way  of  defining  the  average  difference  between  wall 
temperatures  for  uniform  heat  flux.  For  uniform  wall  temperature 
(26)  is  replaced  by 


(to  - tb) 


[(t0-tb)/(<loAoav^X 


av 


which  gives  a Nusselt  number  consistent  with  equation  (25).  For  points 
beyond  the  entrance  region  the  Nusselt  numbers  and  temperature  distri- 
butions are  calculated  as  outlined  in  appendix  C. 


Development  of  Flow  Boundary  Layer 


Mumentum  equations  for  flow  boundary  layers.  - The  momentum  equa- 
tion for  the  fluid  in  a given  space  can  be  written  for  steady  flow  in 


vector  notation  as 


(27) 


or  the  force  acting  on  the  fluid  in  a space  is  equal  to  the  rate  of 
flow  of  momentum  out  of  the  space  (ref.  10,  p.  233).  When  the  X com- 
ponent of  equation  (27)  is  written  for  the  previously  illustrated 
differential  .annulus  (6b  replaced  by  &),  there  results 


(2. 


where  the  differential  of  the  second  integral  is  the  radial  mass  flow 
into  the  annulus  at  6.  The  force  dFx  is  composed  of  pressure  forces 

on  planes  1 and  2 and  a shear  force  at  the  wall.  There  is  no  shear 
force  at  6 because  the  velocity  gradient  is,  by  definition  of  the 
boundary  layer,  zero  at  the  edge  of  the  boundary  layer.  The  force  dFx 


is  then  equal  to 


if(r0-6) 


2J 


dp  - 2rt  Tq  tq  dX 


Substituting  for  dFx  gives,  for  equation  (28), 
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5(CrQ-b)dp  + 


dX  - 2ug, 


Pu(rQ-y)dy 


Pud (rr -y)dy 


’] 


(29) 


In  the  region  outside  the  boundary  layer,  the  flow  is  frictionless,  so 
that 


dp  = - Pu6  dug. 


(30) 


Sub  at  i tut  ing  equation  (30)  in  equation  (29),  rearranging,  and  integrating 
between  the  tube  entrance,  and  X give 


where  the  boundary- layer  thickness  has  been  taken  as  zero  at  X = 0. 
Equation  (31)  is  the  desired  equation  for  flow  in  tubes  relating  the 
boundary -layer  thickness  to  distance  along  the  tube. 


Equations  corresponding  to  equations  (29)  and  (31)  can  be  derived 
for  flow  between  flat  plates.  Those  equations  are 


(33) 
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Equation  (32)  was  first  obtained  by  von  Karman  and  has  been  used  exten- 
sively in  the  incompressible  form  for  calculating  the  growth  of  the 
boundary  layer  for  fluids  flowing  over  surfaces  (ref.  7,  p.  66). 

The  value  of  5 must  be  known  in  terms  of  ug  to  solve  equa- 
tions (31)  and  (33).  This  relation  will  be  obtained  in  the  section 
Conservation  of  mass,  constant  properties. 


Dimensionless  form  of  momentum  equations  for  constant  properties.  - 
Equation  (ol)  for  tubes  can  be  written  in  dimensionless  form  for  constant 
fluid  properties  as 


(34) 


Equation  (34)  can  be  verified  by  substituting  the  definitions  of  the 
dimensionless  quantities  into  the  equation.  Equation  (34)  can  be  put 
into  a form  more  convenient  for  computation  by  using  the  definition  for 
the  differential  of  a product.  Thus 


r0 


d(r0+  %+) 


(W) 


dy+ 


(35) 


Substitution  of  equation  (35)  into  the  last  term  of  equation  (34) 
(the  substitution  can  be  more  easily  made  if  eq.  (34)  is  first  differen- 
tiated with  respect  to  X)  and  simplification  give 
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Equation  (36)  ralataa  X/D  to  b + for  various  values  of  Be 

for  tubes.  The  relation  between  u+  and  y+  is  obtained  from  the 

curve  for  p - 0 in  figure  2.  The  value  of  rQ+  for  a given  Ee  and 

&+  will  be  obtained  in  the  nort  section.  The  corresponding  equation 
for  flow  between  flat  plates  is  found  from  equation  (33)  to  be 


Conservation  of  mass,  constant  properties.  - In  order  to  solve 
equations  (36)  and  (37),  u^+  must  be  known  in  terms  of  Ee.  This 

relation  can  be  obtained  from  the  law  of  conservation  of  mass.  For  a 
constant-area  duct  and  constant  density,  the  law  of  conservation  of  mass 
gives 


UI  = ub 


(38) 


or  the  bulk  velocity  is  constant  and  equal  to  the  uniform  initial 
velocity.  For  a tube,  equation  (38)  can  be  written 


ui 


- y)dy 


u(rQ  - y)dy  + — g-  (r  - 6)2 
r0 


6 


(39) 
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Equation  (30)  can  be  written  in  dimension  lee  a form  as 


y+)dy4 


2u6+(r0+  - 6+) 


(40) 


Tho  equation  corresponding  to  oquation  (40)  for  flow  between  flat  plates 
la 


Re 


+ 4ub+(r0+ 


b+) 


(41) 


Equation  (10)  or  (41)  gives  tho  rotation  between  Re  .and  &+  or  ub+ 
for  a given  value  of  r,,4 *. 


Nuaselt  numbers,  Reynolds  numbers,  and  friction  factors;  constant 
fluid  properties.  - Nusselt  numbers  for  a developing  velocity  distribu- 
tion as  well  as  a developing  temperature  distribution  can  be  calculated 
from  equations  (10)  .and  (12)  or  (14)  and  (16)  as  before.  If  the  assump 
tion  is  made  that  the  thermal-  and  f low-boundary-layer  thicknesses  are 
equal 

and  for  constant  fluid  properties 


(&h+  a b+  for  Pr  - 1 ) ‘! , equation  (12)  for  tubes  can  be  simplified, 


u + (r, 


(j3  = 0)  becomes 
■+)dy+  + - 


V 


(rO+  - &h+) 


'-b 


(42) 


> + - y+)dy+ 


(p0+  - 6h+)‘ 


A similar  equation  cun  bo  obtained  for  flow  between  flat  plates. 


Reynolds  numbers  for  flow  in  tubes  and  between  plates  can  be  calcu- 
lated from  equations  (40)  and  (4l).  For  constant  fluid  properties  the 
Reynolds  number  is  constant  along  the  length  of  the  tube. 


Tho  pressure  drop  between  the  entrance  and  a given  section,  which  is 
required  to  calculate  the  friction  factor,  is  obtained  by  integrating 
equation  (b0). 


Pi 


P u f 


(43) 


Or,  in  dimensionless  form, 

P*  - 


(44) 


4This  assumption  is  used  only  for  determining  tb+  for  the  case  where 

the  Initial  velocity  and  temperature  distributions  are  uniform. 
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for  tubes,  or 

P'  = 8(u6+  r0T  - \ Be6 

(45) 

for  flat  plates . 

The  pressure-drop  parameter  p' 

in  these  equations  is 

defined  as 

(Pi  - P)PD2 
P'  --  o 

(46) 

P2 

The  average  friction  factor  baaed  on  the  pressure  drop  for  a given 
is  defined  by 


f 


av 


_ D(Pj  ~ P) 
2XPub2 


X/D 

(47) 


or 


‘av 


( X/D) Re 


(48) 


This  equation  includes  momentum  changes  caused  by  velocity  profile 
development . 


The  local  friction  factor  based  on  the  local  pressure  gradient  is 
defined  by 


f = . D[d£/d Q 

1 O 

2PUb" 


(49) 


_ 1 dp '/d (X/D) 
2 Re*" 


(50) 


When  f is  to  be  found  from  equation  (50),  the  curve  for  p‘  against  x/D 
can  be  differentiated;  or  the  curves,  calculated  from  equation  (36)  or  (37), 
for  U5+  rQ+  against  x/D  can  be  differentiated  and  dp'/d(X/D)  cal- 
culated from  equation  (44).  The  latter  procedure  proved  to  be  the  more 
accurate . 


The  local  friction  factor  based  on  the  local  shear  stress  at  the 
wall  and  excluding  momentum  changes  due  to  a developing  velocity  profile 
is  defined  by 


2T0 
P ub2 


(51) 
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or 


f-t  - ~^~2  (52) 

V 

and,  since  Re  is  constant  along  the  tube  and  equal  to  2u-j3+  Tq+, 

8r(j+2 

t\  = —V  (53) 

Re1' 

where  the  relation  between  r0+  and  Re  for  a given  6 + is  obtained 
from  equation  (4u)  or  (4l). 

RESULTS  AND  DISCUSSION 

Heat  Transfer  with  Uniform  Heat  Flux,  Uniform  Initial 
Temperature  Distribution,  Fully  Developed  Velocity 
Distribution,  and  Variable  Fluid  Properties 
f ' r Gases;  Frandtl  Number  of  0.73 

Development  of  thermal  boundary  layer.  - Figure  3 illustrates  the 
♦ development  of  the  temperature  profile  in  a tube  for  uniform  heat  flux 

and  a fully  developed  velocity  profile.  The  distance  along  the  tube 
for  various  thermal-boundary- layer  thicknesses  is  calculated  from  equa- 
tion (4)  and  the  temperature  distribution  in  figure  1.  The  temperature 
distribution  in  the  thermal  boundary  layer  Is  obtained  from  figure  1; 
the  temperature  outside  the  boundary  layer  is  uniform,  inasmuch  as  no 
heat  penetrates  that  region.  The  Reynolds  numbers  are  calculated  from 
equations  (ll)  and  (13).  Figure  3 indicates  that  for  a fixed  tube 
diameter  at  a given  distance  from  the  entrance,  the  boundary- layer 
thickness  is  less  for  high  Reynolds  numbers  than  for  low  ones,  or  the 
distance  required  for  a given  degree  of  temperature  profile  development 
is  greater  for  high  Reynolds  numbers  than  for  low  ones. 

Although  the  temperature  distributions  in  figure  3 are  for  0 = 0, 
or  constant  fluid  properties,  similar  distributions  are  obtained  for 
variable  fluid  properties  with  largo  amounts  of  heat  transfer,  as  may 
be  inferred  from  figure  4.  The  variation  of  the  thermal  entrance 
length,  defined  as  the  length  required  for  the  thermal  boundary  layers 
from  the  opposite  walls  to  meet  at  the  center  of  the  passage,  with 
Reynolds  number  and  heat  flux  qi ' is  shown  in  figure  4.  The  quantity 

q^ ' was  calculated  from  equation  (22).  The  quantity  q^ ' rather  than 

0 is  used  as  a heat -flux  parameter  inasmuch  as  it  is  constant  along 
the  length  of  the  tube.  For  the  range  of  qj  ' shown,  t0/tt  varies 

from  about  0.4  to  3.  Figure  4 indicates  that  heat  flux  haa  but  a 
slight  effect  on  the  thermal  entrance  length. 
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Local  Nusselt  numbers  for  tuber,.  - Local  Nusselt  numbers  hD/k-^ 
with  properties  evaluated  at  the  inlet  temperature  are  plotted  against 
X/D  for  various  values  of  Re^  and  q^ ' (heat  flux)  in  figure  5.  The 

heat-tranafor  coefficient  in  the  Nusselt  number  is  based  on  the  differ- 
ence between  wall  and  bulk  temperature  (see  appendix  A).  For 
obtaining  the  curves  Nusselt  numbero  wore  calculated  from  equations  (10), 
(2d),  and  (21);  Reynolds  numbers,  from  equations  (ll),  (18),  and  (19); 
q 1 , from  equation  (22);  and  X/D , from  equation  (4)  for  various  values 

_.f  the  parameters  and  p.  Interpolation  was  necessary  for 

obtaining  the  variation  of  Nu^  with  X/D  for  constant  values  of  Re^ 
and  q^_  ' . 

The  curves  In  figure  5 indicate,  as  might  bo  expected,  that  the 
Nuaselt  numbers  (bet -transfer  coefficients)  close  to  the  entrance  are 
very  high  in  Comparison  with  the  fully  developed  values  farther  down  the 
tube  at  the  ends  of  the  curves.  The  'w  values,  are  high  because 
the  thermal  boundary  layer  is  thin  and  the  temperature  gradients  conse- 
quently are  severe  near  the  entrance.  At  X/D  = 0 the  boundary- layer 
thickness  t>h  is  zero,  so  that  the  heat-transfer  coefficient  is  infinite 

at  the  entrance.  For  uniform  heat  flux  this  means  that  the  temperature 
difference  tu  - tti  must  be  zero  at  X/D  --  0;  for  a finite  temperature 

difference  at  X/D  = 0 the  heat  transfer  per  unit  area  could  not  be 
uniform,  since  it  would  be  infinite  at  the  entrance. 

Inspect  L.  n ->£  the  curves  in  figure  5 chows  that  the  Nusselt  numbers 
very  nearly  r-'U'di  Ui>  ir  fully  developed  values  long  before  the  thermal 
1.  unlasy  layers  h UVe  met  at  the  ends  of  the  curves.  For  instance,  for 

0 (constant  fluid  p report i op)  and  Re^  ~ 100,000,  the  Nusselt  number 

at  X/l‘  -f  9 is  within  2 percent  oi  its  fully  developed  value  at  X/D 
of  lw.-..  This  indicates  that,  for  practical  purposes,  the  thermal  entrance 
lengths  can  %•  considered  to  be  less  than  half  the  values  given  in 
.figure  ' . 


Comparison  -f  the  curves  for  qi ’ =0,  G.0u4,  and  -u.0025  indicates 
a docreano  in  Nuj.  at  a glvon  X/D  for  heat  addition  to  the  gas  and 

an  increase  for  heat  extraction  from  the  gas.  For  both  heat  addition 
and  extraction,  Nuj_  increases  with  X/D  for  largo  values  of  X/D. 

Those  brands  arc  clearly  shown  in  figure  6,  where  Nup  is  plotted  for 
q.|  ’ = o.oO'l  and  -0.0025  for  large  values  of  X/D.  For  q^  ’ =■  0 (constant 

properties)  the  Nusselt  numbers  are  constant  beyond  the  entrance  length. 
The  Nusselt  numbers  for  values  of  X/D  boyond  the  entrance  region  were 
calculated  by  the  method  described  in  appendix  C.  For  heat  extraction 
from  the  gas  the  curvoo  are  cut  at  the  points  shown  because  the  wall 
temperature  roaches  absolute  zero  at  those  points.  The  portions  of  the 
curvos  for  which  the  temperature  is  below  the  liquefaction  temperature 
of  the  gas  should,  of  course,  be  disregarded  because  the  properties  no 
longer  have  the  assumed  variations  with  temperature.  The  increases  in 
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Nusselt  number  with  x/D  are  caused  by  the  variation  of  properties  along 
the  tube  and  across  the  tube.  The  increases  are  due  principally  to  the 
variation  of  t0/tb,  which  decreases  along  the  tube  for  both  heat  addition 

to  and  extraction  from  the  gas. 

It  is  clear  from  figure  6 that  fully  developed  heat  transfer,  in 
the  sense  that  the  heat-transfer  coefficient  becomes  independent  of  X/D, 
cannot  be  obtained  when  the  properties  are  variable.  Fully  developed  heat 
transfer  can,  however,  be  obtained  in  the  sense  that  the  relation  between 
the  Nusselt  number  and  Reynolds  number  becomes  independent  of  X/D  for 
sufficiently  large  values  of  X/D  when  the  properties  are  evaluated  at 
the  proper  reference  temperature,  as  shown  in  figure  7.  The  curves  in 
figure  7 were  obtained  by  cross-plotting  the  curves  for  q_^  1 = 0 from 

figure  4.  They  can  be  used  for  q.^  ’ ■£  0 by  evaluating  the  properties, 

including  density,  in  the  Nusselt  number  and  Reynolds  number  at  the 
reference  temperatures  tx  = x(tQ  - t^,)  + given  in  figure  8.  The 

values  of  x were  computed  by  selecting  values  of  tx  such  that  the 
values  of  Nux  and  Rex,  calculated  from  the  curves  in  figure  5,  from 
equations  (18)  to  (Cl_),  and  from  the  definition  of  tx,  fall  on  the  lines 
for  q.^ ' = 0 in  figure  7.  The  curves  in  figure  8 indicate  that  the 

reference  temperature  for  points  close  to  the  entrance  is  slightly 
closer  to  the  bulk  temperature  than  the  reference  temperature  for  fully 
developed  heat  transfer,  although  its  variation  with  X/D  is  not  large. 

Effect  of  variation  of  heat  transfer  across  boundary  layer  on 
Nusselt  number  for  small  and  constant  fluid  properties.  - For  large 

values  of  bh+  for  fully  developed  flow),  the  effect  of  the  radial 

variation  of  heat  transfer  on  temperature  distributions  was  checked  in 
reference  2 and  found  to  be  small.  However,  for  low  values  of 
such  as  occur  very  close  to  the  entrance,  the  flow  is  nearly  laminar, 
so  that  the  temperature  profile  in  the  boundary  layer  is  not  flat  as  it 
is  for  higher  values  of  b^+.  The  temperature  distribution  in  the 

thermal  boundary  layer  for  laminar  heat  transfer  from  a flat  plate  can 
be  closely  represented  by  a cubic  parabola  (ref.  7,  p.  89),  which  means 
that  the  heat  transfer  per  unit  area  varies  as  1 - y2/^  > as  shown  in 

appendix  B.  This  variation  of  heat  transfer  is  used  herein  for  turbulent 
heat  transfer  for  small  6^+  in  order  to  obtain  the  effect  of  variation 

of  t/Tq  on  the  Nusselt  numbers  in  the  entrance  region.  The  expression 
for  the  temperature  distribution  with  this  variation  of  q./q.Q  is  obtained 
in  appendix  B. 
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The  variation  of  Nu  with  X/D  for  variable  heat  transfer  per 
unit  area  across  the  boundary  layer  (q/q0  = 1 - y2/&h2)  is  shown  in  fig- 
ure 9.  For  comparison,  the  curves  from  figure  5(a)  for  constant  heat 
transfer  (q/q0  = l)  are  shown  by  dotted  lines.  The  curves  in  figure  9 

indicate  that  the  variation  of  heat  transfer  per  unit  area  across  the 
boundary  layer  has  a negligible  effect  on  the  relation  between  Nusaelt 
number  and  X/D,  so  that  the  use  of  the  temperature  distributions  in 
figure  1 appears  Justified,  even  for  low  values  of 

Wall  and  bulk  temperature  distributions  for  tubes.  - Wall  and 
bulk  temperatures  along  a tube  with  uniform  heat  transfer  as  calculated 
from  equations  (23)  and  (24)  ure  plotted  in  dimensionless  form  in 
figure  lo.  As  mentioned  previously,  the  wall  and  bulk  temperatures  are 
equal  at  the  entrance  for  uniform  heat  transfer.  Tho  bulk  temperature 
varies  linearly  with  X/D  inasmuch  as  the  specific  heat  is  assumed 
constant.  For  q^ ' = C the  wall  and  bulk  temperature  curves  are 

parallel  beyond  the  entrance  region.  For  both  heat  addition 

(q^ ' = 0.U04)  and  heat  extraction  (q^ ' = --.0025)  from  the  gas,  the 

wall  and  bulk  temperature  curves  converge  beyond  the  entrance  region 
because  of  the  increase  in  heat-transfer  coefficient  with  X/D  for 
those  portions  of  the  curves. 

Average  Nusaelt  numbers  for  tubes.  - The  variation  of  average 
Nusselt  number  Nuj  av  with  x/D  as  calculated  from  equation  (25)  is 

shown  in  figure  11.  Tho  trends  are  similar  to  those  for  local  Nusselt 
numbers  but  the  changes  are  more  gradual. 


The  relation  between  Nux  ,iV  and  Rex  ay  for  various  values  of 

X/D  is  given  in  figure  12.  The  curves  were  obtained  by  cross-plotting 
the  curves  for.  q^ ' - 0 in  figure  11 (a)  and  can  be  used  for  q^ ' / 0 

by  evaluating  the  properties  in  the  Nusselt  number  and  Reynolds  number 
at  tx,av  " >:(t0,av  " H.av^  + H, av>  vheru  x is  given  in  figure  13. 

The  values  of  x wore  computed  by  selecting  values  of  tx  ay  such  that 
the  values  of  Nux  av  and  Rex  av  calculated  as  outlined  in  appendix  C 
fall  on  the  lines  for  q^ ' = 0 in  figure  12.  Figure  13  indicates  that 

the  value  of  x increases  continuously  with  X/D  for  both  heat  addition 
to  and  extraction  from  the  gas.  These  trends  are  indicated  experimentally 
for  heat  addition  in  reference  11  . whore  a constant  reference  temperature 
used  for  a range  of  values  of  X/D  was  found  to  give  average  Nusselt 
numbers  for  high  heat  fluxes  at  low  values  of  X/D  which  were  slightly 
higher  than  those  for  low  heat  fluxes.  The  values  of  x given  in 
figure  13  indicate  that  for  large  values  of  X/D,  x>1,  so  that  for  heat 
addition  tXjav  is  higher  than  the  average  wall  temperature  although 

not  necessarily  higher  than  the  exit  wall  temperature.  The  variation 

of  x with  X/D  is  caused  by  the  variation  of  properties  along  the 
tube . A more  satisfactory  correlation  would  perhaps  be  obtained  by 
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using  x only  to  account  for  the  variation  of  properties  across  the 
tube  and  using  another  factor  to  account  for  the  variation  along  the 
tube.  The  curves  in  figures  12 'and  13  might,  however,  provide  a con- 
venient means  of  utilising  the  results  of  the  analysis  for  some  purposes. 

Nusselt  numbors  and  temperature  distributions  for  flow  between  flat 
plates . - Nusselt  numbers,  wall  temperature  distributions,  and  bulk  tem- 
perature distributions  for  flow  between  flat  plates  with  uniform  heat 
transfer  and  fully  developed  velocity  distribution  are  plotted  in 
figures  14  to  16.  The  curves  agree  closely  with  the  corresponding  curves 
for  flow  in  a tube.  For  Instance,  the  curves  for  local  Nusselt  numbers 
plotted  in  figure  14  are  between  3^  and  5 percent  higher  in  the  fully 

developed  portions  of  the  curves  than  the  corresponding  curves  for  tubes 
in  figure  5.  The  portions  of  the  curves  close  to  the  entrance  agree 
oven  more  closely.  This  agreement  indicates  that  results  for  turbulent 
heat  transfer  in  tubes  can  be  used,  with  small  error,  for  flat  plates 
when  an  equivalent  diameter  equal  to  twice  the  plate  spacing  is  used. 

It  should  be  emphasized  that  in  applying  the  results  of  the  fore- 
going analysis  the  same  assumptions  regarding  the  variation  of  properties 
with  temperature  should  be  made  as  were  made  in  the  analysis,  that  is, 
constant  specific  heat,  and  viscosity  and  thermal  conductivity  both 
proportional  to  tu-68. 

It  should  alno  be  mentioned  that  the  results  in  this  analysis  do 
not  include  the  low  Pec  let  number  effect  described  in  reference  8 inas- 
much as  the  temperature  distributions  in  figure  1 do  not  include  that 
effect.  The  Nusselt  numbers  at  Rej^  = 10, GOO  may  therefore  be  slightly 

high,  but  the  ratios  of  local  to  fully  developed  Nusselt  numbers  should 
be  accurate . 

Effect  of  frictional  heating  on  entrance  length.  - The  effect  of 
frictional  heating' on  the  x/d  for  a given  can  be  approximately 

accounted  for  by  replacing  t^+  and  t+  by  T$+  and  T+,  respectively, 
In  equation  (4)  and  setting  p/pQ  = l/(l  - (3T+  - au+2)  (ref.  9).  In 

using  that  procedure  It  is  assumed  that  the  total-temperature  gradient  is 
zero  at  the  edge  of  the  thermal  boundary  layer,  a good  assumption  for 
Prandtl  numbers  close  to  unity. 

The  variation  of  thermal  entrance  length  (x/D)^e  with  Reynolds 

number  and  heat  flux  with  frictional  heating  (a  = 0.00025)  is  shown 
in  figure  17.  Mach  numbers  for  air  at  the  tube  center  at  (x/D)he  are 

indicated  in  the  figure.  The  Mach  numbers  were  calculated  from 
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which  can  bo  obtained  from  the  definitions  of  M,  Tc,  uc+,  and  a.  Com- 
parison of  figure  17  with  figure  4 indicates  that  frictional  heating 
increases  the  entrance  length,  although  the  increase  is  not  large  for 
subsonic  flow.  Inasmuch  as  Nusselt  numbers  for  fully  developed  heat 
transfer  are  not  affected  by  frictional  heating  (ref.  9),  the  effect 
of  frictional  heating  on  the  relations  between  Nusselt  number  and  X/D 
probably  is  slight  when  the  total  rather  than  the  static  bulk  tempera- 
ture is  used  in  the  definition  of  the  heat-transfer  coefficient. 


Heat  Transfer  with  Uniform  Wall  Temperature,  Uniform  Initial 
Temperature  Distribution,  Fully  Developed  Velocity 
Distribution,  and  Constant  Fluid  Properties  for 
Gases;  Prandtl  Number  of  0.73 

Relations  among  Nussolt  number,  Reynolds  number,  and  X/D  for  the 
case  of  tubes,  as  calculated  from  equations  (9)  to  (11),  are-  plotted  in 
figure  18.  Co-mparis.n  of  figure  18  with  figure  5 indicates  that  the 
Nusselt  numbers  for  uniform  wall  temperature  are  very  slightly  lower 
than  those  for  uniform  heat  flux.  The  heat  transfer  per  unit  area  qsj 

must  be  infinite  at  X/D  - 0 for  uniform  wall  temperature.  This  occurs 
only  for  an  infinitesimal  distance  dX,  however,  so  that  a finite  amount 
of  heat  is  transferred . 


* 
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Heat  Transfer  with  Uniform  Heat  Flux,  Uniform  Initial  Velocity  and 
Temperature  Distributions,  and  Constant  Fluid  Properties 
for  Gases;  Prandtl  Number  of  0.73 

For  calculating  the  variation  of  Nusselt  number  with  X/D  for 
this  case  the  f low-boundary-layer  thickness  as  well  as  the  thermal- 
boundary-layer thickness  should  be  known  for  various  values  of  X/D. 
However,  in  the  calculation  of  Nusselt  numbers,  it  is  considered 
sufficiently  accurate  to  assume  the  thicknesses  of  the  thermal  and 
flow  boundary  layers  to  be  equal,  inasmuch  as  the  Prandtl  number  is 
close  to  unity. 

Relations  among  Nusselt  number,  Reynolds  number,  and  X/D  for 
tubes,  as  calculated  from  equations  (10),  (42),  (40),  and  (4),  are 
plotted  in  figure  19.  Comparison  of  figure  19  with  figure  5 indicates 
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that  the  Nusselt  numbers  for  a uniform  Initial  velocity  distribution  are 
higher  than  those  for  a fully  developed  velocity  distribution.  This 
result  might  be  expected  because  of  the  higher  friction  in  the  case  of 
the  uniform  Initial  velocity  profile. 

Average  Nusselt  numbers  as  calculated  from  equation  (25)  are  plotted 
in  figure  20.  Figure  20(b),  which  is  a logarithmic  plot,  compares  the 
predicted  relations  between  Nuav  and  X/D  with  the  empirical  relation 

Nu  a(x/D)"°-1.  Curves  having  a slope  of  -0.1  on  the  log-log  plot  agree 
closely  with  the  predicted  lines  for  X/D  between  6 and  60.  The  agree- 
ment of  the  data  given  in  reference  11  with  the  empirical  relation  for 
values  of  N/D  as  high  as  120  might  be  due  to  the  use  of  a constant 
rather  than  a variable  reference  temperature  as  given  in  figure  13. 

Nusselt  numbers  for  flow  between  flat  plate3  are  given  in  figure  21. 
Comparison  of  figure  21  with  figure  19  indicates  very  close  agreement 
between  the  Nusselt  numbers  for  tubes  and  flat  plates. 


Heat  Transfer  with  Uniform  Wall  Temperature,  Uniform  Initial 
Velocity  and  Temperature  Distributions,  and  Constant  Fluid 
Properties  for  Gases;  Prandtl  Number  of  0.73; 

« /and  Comparison  with  Experimental  Data 

The  variation  of  Nusselt  number  with  X/D  for  this  case  is  given 
in  figure  22.  The  curves  agree  closely  with  those  for  the  corres- 
ponding case  for  uniform  heat  flux  given  In  figure  19.  The  Nusselt 
numbers  for  uniform  wall  temperature  are  slightly  lower  than  those  for 
uniform  heat  flux. 

A comparison  between  analytical  and  experimental  results  for  uni- 
form wall  temperature  and  uniform  Initial  velocity  and  temperature  dis- 
tributions is  given  in  figure  23.  The  experimental  curves  represent 
data  from  reference  12  for  air  flowing  in  a tube  at  uniform  temperature 
with  a bellmouth  entrance  and  a screen  at  the  entrance  and  should  repre- 
sent the  case  analyzed  in  this  section.  The  screen  at  the  entrance 
should  insure  a turbulent  boundary  layer  throughout  the  tube.  The 
figure  shows  that  the  agreement  betwoen  analytical  and  experimental 
results  for  this  case  is  very  good. 

The  lack  of  a comparison  of  the  present  analytical  results  for  a 
uniform  wall  temperature  and  a fully  developed  velocity  distribution 
with  the  data  from  reference  12,  which  were  taken  with  a long  approach 
section  at  the  entrance,  should  perhaps  be  explained.  It  appears  that 
the  sharp  edge  at  the  beginning  of  the  long  approach  section  produced 
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large  disturbances  which  were  not  damped  out  In  the  length  of  approach 
section  used.  The  data  in  reference  12  for  heat-transfer  coefficients 
with  a sharp  entrance  without  an  approach  section  Indicate  such  results. 
For  that  case  the  heat-transfer  coefficients  are  higher  throughout  the 
tube  than  thoso  obtained  with  a bellmouth  entrance. 

Average  Nusselt  numbers  calculated  from  equation  (25)  are  plotted 
in  figure  24.  As  in  the  case  of  local  Nusselt  numbers,  the  average 
Nusselt  numbers  for  a uniform  wall  temperature  in  figure  24  are  slightly 
lower  than  those  for  a uniform  hoat  flux  in  figure  20. 


Heat  Transfer  with  Uniform  Heat  Flux,  Uniform  Initial 
Temperature  Distribution,  Fully  Developed  Velocity 
Distribution,  and  Constant  Fluid  Properties  for 
Liquid  Metals;  Prandtl  Number  of  0.01 

Temperature  distributions  for  fully  developed  flow  calculated  as 
described  in  reference  8 and  plotted  in  figure  25  were  used  in  the 
thermal  boundary  layer  for  calculating  the  variation  of  Nusselt  number 
with  X/D  for  liquid  metals.  The  velocity  distributions  were  obtained 
frem  the  curve  for  p = o in  figure  2.  The  curves  for  Nu  against 
X/D  in  figure  26  were  calculated  from  equations  (4)  and  (10)  to  (12)  * 

(p  = o) . An  unexpected  feature  of  these  curves  is  the  slight  increase 
of  Nusselt  number  with  X/D  at  large  values  of  X/D  and  Reynolds  num- 
ber. This  increase  is  due  to  the  change  in  shape  of  the  temperature 
profile  with  &k+  for  low  Prandtl  numbers.  Thus,  as  increases, 

the  temperature  gradient  at  the  wall  increases  slightly  in  some 
instances  because  of  the  added  turbulence  in  the  boundary  layer  at 
higher  values  of  6^+.  This  effect  is  indicated  experimentally  by 

the  data  in  reference  13. 

A plot  of  Nusselt  number  against  Pe/(X/D)  for  various  Peclet 
numbers  is  given  in  figure  27.  If  the  heat  transfer  contributed  by 
eddy  diffusion  were  neglected,  the  curves  for  various  Reynolds  numbers 
would  fall  essentially  on  a single  line  as  in  reference  5.  Examination 
of  the  curves  in  figure  27  indicates  that  for  Reynolds  numbers  below 
1000  the  curves  could  be  represented  approximately  by  a single  line. 

This  line  is  slightly  higher  than  the  one  obtained  in  reference  5 inas- 
much as  a uniform  wall  temperature  rather  than  a uniform  heat  flux  was 
postulated  in  that  reference. 
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Friction  Factors  with  Uniform  Initial  Velocity  Distribution 
and  Constant  Fluid  Properties 

Throe  friction  factorn -are  discussed  in  this  section:  fT,  based 

on  the  shear  stress  at  the  wall;  f,  based  on  the  static- pressure  gradient; 
and  fa_,  based  on  the  pressure  drop  along  the  tube.  The  friction  factor 

t\  is  calculated  from  equation  (53);  f,  f^um  equation  (50);  and  fa  , 

from  equation  (48).  The  Reynolds  number:,  aid  values  of  X/D  corres- 
ponding to  these  friction  factors  are  can’.uj a from  equations  (41)  and 
(36).  The  velocity  distributions  in  the  fjov  boundary  layer  are 
obtained  from  the  curve  for  0 - 0 in  figure  2 . 

The  friction  factors  for  a tube  based  on  the  shear-  r tress  at  the 
wall  and  on  the  static  pressure  gradient  are  plotted  in  figures  28  and 
C*2 , respectively.  The  friction  factors  based  on  the  pressure  gradient 
are  higher  than  those  based  on  the  shear  stress  in  the  entrance  region 
because  the  former  Include  the  pressure  decrease  caused  by  the  increase 
in  momentum  of  the  fluid  as  the  velocity  profile  develops  as  well  as 
the  pressure  decrease  caused  by  the  shear  stress  at  the  wall.  For 
fully  developed  flow  the  two  friction  factors  are  identical.  Average 
friction  factors  based  on  the  pressure  drop  are  plotted  in  figure  30 
and  arc:,  of  course,  higher  than  the  local  values. 

Tile  computed  values  of  f and  ffiV  are  subject  to  3ome  inaccuracy 

becau.T.  _f  errors  associated  with  measuring  the  slope  of  a curve  (eq. 

(50))  and  because  the  two  torms  in  equation  (44),  which  'ire  subtracted, 
are  nearly  equal  ciose  to  the  entrance.  Thu  results  presented  should, 
however,  be  well  within  the  limits  of  accuracy'  of  experimentally  deter- 
mined friction  factors. 

Friction  factors  for  flow  between  flat  plates  are  plotted  in 
figures  31  to  55  and  are  found  to  be  very  similar  to  those  for  a tube. 


Comparison  of  Local  Nusselt  Numbers  and  Friction  Factors 

A comparison  of  local  Nusselt  numbers  and  friction  factors  for 
various  cases  is  presented  in  figure  34,  where  the  ratio  of  local 
Nusselt  number  or  friction  factor  to  the  fully  developed  value  is 
plotted  against  X/D  for  a Reynolds  number  of  100,000.  The  quantity 
Nu/Nu^  is  smaller  than  f/f^  for  a given  X/D  in  all  cases  except 

that  of  heat  transfer  to  a liquid  metal  (Pr  = 0.01).  The  larger  values 
of  Nu/Nu^  for  a liquid  metal  are  apparently  caused  by  the  fact  that 

the  temperature  profile  at  low  Prandtl  numbers  is  similar  to  that  for 
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laminar  flow.  At  very  high  Reynolds  numbers,  where  the  dip  occurs  in 
the  curves  for  Nusselt  number  plotted  against  X/D,  and  at  very  low 
Reynolds  numbers,  the  values  of  Nu/Nu^  may  be  lower  than  f/f^  for 

liquid  metals  (fig.  2d).  For  gases  the  values  of  Nu/Nu^  are  lower 

than  f/f^  for  a given  value  of  X/D  because  of  the  momentum  pressure 
loss  for  velocity  profile  development  included  in  the  friction  factor 
f;  there  is  no  counterpart  for  this  momentum  pressure  loss  in  heat 
transfer. 

It.  should  be  mentioned  that  the  preceding  conclusions  may  apply 
only  to  the  case  analyzed,  that  is,  to  the  case  where  the  boundary  layer 
is»  turbulent  throughout  the  tube.  The  turbulent  boundary  layer  might  be 
produced  by  a bellmouth  with  a strip  of  sand  paper  at  the  entrance  or 
by  screens.  For  flow  in  which  the  boundary  layer  is  partially  laminar, 
or  for  which  there  are  large  disturbances  at  the  entrance,  such  as  might 
be  caused  by  a right -angle-edge  entrance,  the  conclusions  might  be 
altered . 


SUMMARY  OF  RESULTS 

The  following  results  were  obtained  from  the  analytical  investiga- 
tion of  the  entrance  region  for  turbulent  heat  transfer  and  flow  in 
smooth  tubes  and  between  parallel  flat  plates: 

1.  The  thin  thermal  boundary  layers  and  consequently  severe  tem- 
perature gradients  at  the  wall  near  the  entrance  produced  high  heat 
transfer  near  the  entrance.  Approximately  fully  developed  heat  transfer 
was,  in  general,  attained  in  an  entrance  length  less  than  10  diameters. 

2.  The  effect  of  variable  properties  on  the  local  Nusselt  number 
correlation  for  gases  with  fully  developed  velocity  profiles  was  to 
decrease  the  Nusselt  .number  with  conductivity  based  on  the  inlet  tem- 
perature for  heat  addition  to  the  gas  and  to  Increase  it  for  heat 
extraction  from  the  gas  at  a given  ratio  of  distance  from  the  entrance 
to  tube  diameter  X/D  and  inlet  Reynolds  number.  For  large  values  of 
X/D,  the  Nusselt  number  with  conductivity  based  on  the  inlet  temperature 
increased  continuously  with  X/D  for  both  heat  addition  to  and  extrac- 
tion from  the  gas.  A heat-tran3fer  coefficient  which  is  independent  of 
x/D  is  therefore  not  obtained  for  heat  transfer  with  variable  properties, 
even  in  the  fully  developed  region. 

3.  The  reference  temperature  used  for  evaluating  the  properties 
iii  order  to  eliminate  the  effect  of  heat  flux  on  the  curves  for  local 
Nusselt  number  against  Reynolds  number  at  various  values  of  X/D  varied 
only  slightly  with  X/D.  The  reference  temperature  for  average  Nusselt 
numbers,  however,  increased  continuously  with  X/D  for  heat  addition 
to  the  gas  and  was  higher  than  the  average  wall  temperature  but  lower 
than  the  exit  wall  temperature  at  very  large  values  of  X/D. 


2857 


NACA  TN  3016 


29 


4.  The  variation  of  Nusoolt  number  with  X/D  for  a gas  with  uni- 
form initial  temperature  and  velocity  distributions  and  uniform  heat  flux 
agreed  with  the  experimentally  determined  relation,  Nu«c (x/D)"^-^,  for 
valuer,  of  X/D  between  6 and  6u. 

5.  The  predicted  variation  of  Nusseit  number  with  X/D  for  air 
with  a uniform  wall  1 < >mperaturo  and  uniform  initial  velocity  and  tem- 
perature distributions  agreed  very  well  with  the  experimental  values 
f r - m re f erence  12 . 

0.  The  variation  of  Hue  jo  it  number  with  Graotz  number  for  liquid 
metals  (Prandti  number  ef  .ei)  was  represented  approximately  by  a single 
line  for  Pec  let  numbers  below  Ida  j. 


7.  The  ratio  .f  fricti  n factor  baaed  on  the  static- pres sure 
gradient  in  the  entrance  region  to  the  fully  developed  value  was  greater 
than  the  corresponding  ratio  of  Nurse It  numbers  for  gases  in  the  entrance 
region  because  of  the  m men turn  pressure  loss  associated  with  the  velocity 
prof  i le  d l vo  1 e pmen t. . 

6.  Essentially  the  same  variations  of  Husselt  number  and  friction 
factor  with  X/D  wor<  obtained  for  flow  between  parallel  flat  plates 
as  wore  obtained  for  a tube  when  the  equivalent  diameter  equal  to  twice 
the  plat<o  spacing  was  used  f or  the  plates. 

L'-wir  Flight  Pr -pulsion  Laboratory 

National  Advisory  Ccmmittee  for  Aeronautics 
"love land,  Ohio,  July  21,  lt-53 
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APPENDIX  A 


SYMBOLS 


The  following  symbols  are  used  in  this  report: 

A area,  sq  ft 

Cp  specific  heat,  of  fluid  at  constant  pressure,  Btu/ ( l'> ) (°F) 

cD  specific  host  „.f  fluid  at  constant  pressure  at  wall, 

0 Btu/(lb)(  F) 


D 

d 


h 


j 


r:b 


i:i 


U 


U 


inside  diameter  of  tube,  ft 

exponent;  value  depends  on  variation  of  viscosity  of  fluid 
with  temperature 

acceleration  due  to  gravity,  /i.2  ft/eec" 

local  heat-transfer  coefficient,  q(l/t - t^,  Btu/(sec)(sq  ft)(°F) 

average  heat -transfer  c efficient,  q, . ttV/(tc,  - t^)ay 

mechanical  equivalent  of  heat,  776  ft-lb/Btu 

therms  j.  conductivity  of  fluid,  Btu/(sec)(sq  ft)(°F/ft) 

thermal  conductivity  of  fluid  evaluated  at  t^, 

Btu/ ( sec ) ( sq  ft ) ( °F/f t ) 

thermal  conductivity  of  fluid  evaluated  at  tj_, 

Btu/ ( so  c ) ( sq  ft ) ( °F/ ft ) 

thermal  conductivity  of  fluid  evaluated  at  t^, 

Btu/ (sec) (sq  ft)(uF/ft) 

thermal  conductivity  of  fluid  evaluated  at  tx, 

Btu/ (sec) (sq  ft) (°F/ft) 


thermal  conductivity  of  fluid  evaluated  at  tx  av, 
Btu/ (sec) (sq  ft)(°F/ft) 
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M Mach  number  at  tuba  canter  at  X/D 

n constant 

P static  pressure 

P|  static  pressure  at  inlet,  lb/aq  ft  abo 

q rate  of  heat  transfer  toward,  tubu  center  per  unit  area, 

Btu/(soe) (oq  ft) 

q rate  of  heat  transfer  at  inside  wall  toward  tube  center  per 

unit  area,  Btu/(oec)(sq  ft) 

q,  t a vera  e rate  ..f  heat,  transfer  at  inside  wall  toward  tube 

"HV  center  per  unit  area  for  a given  X/D,  Btu/(sec)(sq  ft) 

r radius,  distance  from  tube  center,  ft 

r.  , inside  tube  radius  or  one-half  distance  between  plates,  ft 

T absolute  t.tn,  temperature,  ^R 

Tc  absolute  t.t.'ii  temperature  at  tube  center,  UR 

Id  total  temperature  outside  thermal  boundary  layer,  °R 

t absolute  static  temperature,  V,R 

t|.;  bull:  static  temperature  of  fluid  at  cross  section  of  tube,  °R 

tb  av  average  bulk  temperature , *R 

t^  inlet  static  temperature  of  fluid,  lR 

tx  reference  temperature  for  local  Nuoselt  and  Reynolds  numbers, 

x(tu  - tb)  + tb,  °R 

t reference  temperature  for  average  Nuoselt  and  Reynolds  numbers, 

y Cl  V » v 

x'\.,av  " ’"bjav^  + ^bjUV 

tj,  temperature  of  fluid  outside  thermal  boundary  layer 
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t-Q  absolute)  wall  temperature,  °R 

t('  >a v average  wall  temperature , °R 

(t/.-tk)  average  difference  between  wall  and  bulk  temperature,  UR 
u time-average  velocity  parallel  to  axis  of  channel,  ft/sec 

Ujj  bulk  velocity  at  cross  section  of  tube,  ft /sec 

Ui  velocity  of  fluid  at  inlet,  ft/m-c 

u^  velocity  outside  flow  boundary  layer,  ft /sec 


y 

6 


&h 


c 

M 


Mj; 

dx,av 


dr, 


P 


velocity  vector,  ft/sec 
distance  from  entrance,  ft 

number  used  Ln  evaluating  arbitrary  temperature  in  tube, 
t:c  "r  tx,av 

distance  fr_r.  wail,  ft 

flow- boundary  layer  thickness,  ft 

thermal -boundary -layer  thicluess,  ft 

coefficient  of  eddy  diffusivity,  sq  ft/soc 
abs  ..lute  viscosity  of  fluid,  (lb)(sec)/sq  ft 
absolute  viscosity  of  fluid  evaluated  at 

absolute  viscosity  of  fluid  evaluated  at 

absolute  viscosity  of  fluid  evaluated  at 

absolute  viscosity  of  fluid  evaluated  at 
(lb)(sec)/sq  ft 

absolute  viscosity  of  fluid  evaluated  at 
mass  density  of  fluid,  (lb) (sec^O/ft^ 


tb,  ( lb) (sec)/sq  ft 
ti}  (lb) (sec)/sq  ft 
tx,  (lb) (see)/rq  ft 
^xg-iv* 

tQ,  (lb) (sec)/ sq  ft 
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ro 

Cn 

cn 


pb,av 

pi 

p* 

px,av 

P0 


T 


• L' 


mass  density  of 
mass  density  of 
mass  density  of 
mass  density  of 
mass  density  of 
mass  density  of 

shear  stress  in 
shear  stress  in 


fluid  evaluated  at 
fluid  evaluated  at 
fluid  evaluated  at 
fluid  evaluated  at 
fluid  evaluated  at 
fluid  evaluated  at 


tb,  (lb)(sec2)/ft4 

(lb)  (sec2) /ft4 
t^,  (lb) (sec2 )/ft4 
tx,  (lb) (see2)/ft4 
tXjav»  (lb) (sec2)/ft4 
t0,  (lb)  (sec2)/ft4 


fluid,  lb/sq  ft 
fluid  at  wall,  lb/sq  ft 


Dimensionless  .groups: 
f friction  factor 

fav  friction  factor 

fully  developed 
fx  friction  factor 

Gz  Graetz  number,  Re  Pr/(X/D) 

Nu  Nusselt  number,  h D/k 

Nu„,r  average  Nusselt  number  for  constant  properties  defined  by 

a.\  , , 

eq.  (2b) 

Nub  IJ us 30 It  number  with  thermal  conductivity  evaluated  at  tb , 

h D/kb 

Nu^  fully  developed  Nusselt  number 

N^  Nusselt  number  with  conductivity  evaluated  at  t.^ 


based  on  static-pressure  gradient,  -P 

2p  ub"- 
D(p1  - p) 

based  cn  static-pressure  drop,  — 

2 X P ub2 

friction  factor 

A 

based  on  shear  stress  at  wall,  2tq/(p  ubu) 
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Nu^  av  average  Nuaaelt  number  defined  by  eq.  (25)  with  conduc- 
tivity evaluated  at  tj_ 

Nux  Nuaaelt  number  with  conductivity  evaluated  at  tx 

Nux  ay  average  Nuaaelt  number  with  conductivity  evaluated  at  tx  av 

Nu0  Nuaaelt  number  with  thermal  conductivity  evaluated  at  t,j 

p'  pre j3uro-drcp  parameter,  (p^  - p)pD'1/p^ 


Pe 

Pr 

Pr0 

^i  ’ 

Ee 

Beb 

Re^ 

Eex 

Hex,av 

Ee0 


Pec  let  number.  Re  Pr 
Prandtl  number,  Cp  g p/i: 

Prandtl  number  with  properties  evaluated  at  t,j 
boat -transfer  parameter,  qrj/(Cp  g piuiti) 

Reynolds  number,  p ub  D/p 

Reynolds  number  with  density  and  viscoBity  evaluated  at  tv, 

v*  ^ 

Reynolds  number  with  density  and  viscosity  evaluated  at  t^, 
piubD 
Pi 

Reynolds  number  with  density  and  viscosity  evaluated  at  tx, 
Pxub  D/px 

Reynolds  number  with  density  and  viscosity  evaluated  at  tx  av, 

„ px,av  pi 

Re<  — J 

pb,av  Px,av 

Reynolds  number  with  density  and  viscosity  evaluated  at  tp, 
p0ub  D/p0 

tube  radius  parameter,  7 r,v 

P0/p0  0 
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ro 

CD 

oi 


total-temperature  parameter, 


(tQ  - T)c  g T0  1 - T/t( 


ao"^To/  po 


Mi  -- 

p \ t 


temperature  parameter , 


(t0  - t)c  g 1 - t/t0 


^oV^o/pc 


bulk-temperature  parameter, 


1 1 . ^ 
P V t0 


1 - T” 


velocity  parameter, 


u 


bulk -velocity  parameter, 


ub 


^Tq 


value  of  u+  at  tube  center 
u. 


-y+ 

a 

P 

5 + 


he 


thermal  entrance  length  divided  by  diameter 


Wi/Pfi 

wall  distance  parameter,  ^ — h y 

P0/p0 

compressibility  or  frictional  heating  parameter,  '*0/(2cJcptopO^ 
heat-transfer  parameter,  <1q Vto/p0  / (e- 


t 


0 0 } 


dimensionless  flow-boundary- layer  thickness. 


dimensionless  thermal -boundary- layer  thickness, 


p0 

■Wp0 


An  integral  sign  having  open  brackets  for  the  upper  limit 

indicates  that  the  variable  of  integration  should  be  used  as  the 
upper  limit. 
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APPENDIX  B 


ANALYSIS  INCLUDING  EFFECT  OF  VARIATION  OF  HEAT  TRANSFER 
ACROSS  BOUNDARY  LAYER  FOR  SMALL  bh+  WITH 
CONSTANT  FLUID  PROPERTIES 


As  was  mentioned  in  the  section  Effect  of  variation  of  heat  transfer 
across  boundary  layer  on  Nusselt  number  for  small  b^*  and  constant 
fluid  properties,  the  effect  of  variation  of  tAq  across  the  boundary 
layer  for  small  bj1+  should  be  checked.  For  laminar  flow  over  a flat 
plate  the  temperature  distribution  in  the  thermal  boundary  layer  can  be 
closely  approximated  by  a cubic  parabola.  For  a cubic  parabola, 
l/lQ  = 1 - //bjj2  as  shown  in  the  following: 


For  laminar  flow, 

4 = - >'  dt/dy  = q0(l  - y2/^2)  (Cl) 


Integration  of  equation  (Cl)  gives 


y 


- Xg  - £ (to  - t) 

V 


36h2  To 


(C2) 


Evaluation  of  equation  (C2)  at  the  edge  of  the  boundary  layer  gives 

^ (t0  - t6)  =|bh  (C3) 

Combination  of  equations  (C2)  and  (C3)  results  in 


(C4) 


which  is  the  same  as  the  equation  on  page  89  in  reference  7.  The 
assumption  that  q/qQ  = 3 - y2/®h2  therefore  leads  to  a cubic  parabola 
for  the  temperature  distribution  for  laminar  flow.  The  same  assumption 
is  used  herein  for  turbulent  flow  for  small  b^"1".  The  equation  for 
turbulent  heat  transfer  can  be  written 
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q - - (k  + cp  g pe)  dt/dy  (C5) 

With  the  assumption  for  q/qg  and  c = n2uy  for  flow  close  to  the  wall 
(ref.  6),  equation  (C5)  becomes,  in  dimensionless  form. 


+2 


■ [b  * 


dt+ 

dy+ 


(C6) 


where  t+ ' represents  the  temperature  distribution  with 
2/2 

q/q^  = 1 - y+  /b^+  as  distinguished  from  t+,  which  is  the  dimensionless 
temperature  parameter  for  q/q^  = 1.  Separation  of  variables  and 
integration  of  equation  (C6)  result  in 


t+' 


Iz. 


1 2 + + 
P?  + n u V 


Iz. 


n2u+y+ 


(C7) 


The  first  term  on  the  right  side  of  equation  (C7)  is  t+  for  constant 
properties  (ref.  6),  ao  that 


a t+ . _i_  ry+  7+z  *7+ . . 

+2  | JL  + n2u+y+ 

h Jo 


(ce) 


The  values  of  t+  are  obtained  from  figure  1.  For  obtaining  the 
Nusselt  numbers  the  values  of  t+ ' from  equation  (C8)  are  used  in 

equations  ( 12 ) and  (4)  in  place  of  t+ . 
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APPENDIX  C 


CALCULATION  OF  NUSSELT  NUMBERS,  REYNOLDS  NUMBERS,  AND  TEMPERATURE 
DISTRIBUTIONS  FOR  FLOW  BEYOND  ENTRANCE  REGION 


For  points  beyond  the  entrance  region  the  relation  between  Nusselt 
number  and  Reynolds  number  with  the  properties,  Including  density, 
evaluated  at  the  reference  temperature  tx  - x(tQ  - t^)  + t-^,  where 

x = 0.4  for  heat  addition  and  0.6  for  heat  extraction  from  the  gas, 
can  be  obtained  from  the  curve  for  fully  developed  flow  In  figure  7. 
From  equations  (18)  to  (21)  and  the  definition  of  tx  it  can  be  shown 

that 


and 


Nui 


(Bl) 


(B2) 


where  d = 0.68  for  air  and  most  other  gases.  The  value  of  Nuj  for 
a given  value  of  Re^  can  be  obtained  from  equations  (Bl)  and  (B2) 
provided  the  values  of  t^/t^  and  t0/t^  are  known.  The  quantities 
tjj/t^  and  t^/tjj  can  be  written  in  terms  of  the  quantities  plotted  in 
figure  10  as 


(Bo) 


(see  eq.  (23))  and 


t0  v 

(X1) 

i 

+ i 

tb 

r/V  ti\ 

i 

V 

[In ) 

^i’ 

+ i 

(B4) 


In  the  procedure  for  obtaining  the  value  of  the  quantity  in 

brackets  in  the  numerator  is  first  determined  by  extrapolating  a known 
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curve  for  that  quantity  plotted  against  x/D.  The  value  of  Nu^  can 
then  be  obtained  from  equation  (B2).  With  that  value  of  Nuj_  a nev  and 


more  accurate  value  of  [(VMAi] 
tion  (24) . The  calculation  of  Nu^ 


— can  be  calculated  from  equa- 
^i 

can  then  be  repeated  using  the  new 


value  of 


[(t0-t1)/ti] 


Average  Nusselt  numbers  with  the  conductivity  evaluated  at  the 
inlet  temperature  Nu^  av  can  be  calculated  from  equation  (25),  where 

the  values  of  Nu^  are  calculated  as  outlined  previously.  It  is 

desirable  to  calculate  average  Nusselt  numbers  and  Reynolds  numbers  with 
properties  evaluated  at  the  average  bulk  or  the  average  wall  temperature 
for  given  values  of  x/D.  The  average  bulk  temperature  rise  for  a given 
X/D  is  one-half  the  temperature  rise  at  X/D,  since  the  specific  heat 
is  assumed  constant.  Therefore, 


or 


(B5) 


From  the  definitions  of  hay,  Nuj_  av,  Re1,  and  Pr,  it  can  be  shown  that 


^0,av  ^b,av  ^i  Be^  Pr 

- - * -f  ■ ’ 

.1.  -L  7 _ 


Nu 


i,av 


(B6) 


Average  Nusselt  numbers  and  Reynolds  numbers  with  properties  evaluated 
at  ^b  av>  ^av'  or  art)itrary  temperature 

^x,av  = x^0,av  “ ^b,av^  + ^b,av 
can  be  calculated  by  use  of  equations  (B5)  and  (B6). 
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Figure  5.  - Concluded.  Variation  of  Nusselt  number  IlUj  with  x/D  and  Reynolds  number  for  gas  flowing  in 
a tube.  Pr,  0.73.  Uniform  heat  flux,  uniform  initial  temperature  distribution,  and  fully  developed 
velocity  distribution. 


Figure  6.  - Concluded.  Variation  of  Nusselt  number  Nu^  with  X/D  and  Reynolds 
number  for  gas  flowing  in  a tube  for  large  values  of  X/D..  Pr,  0.73.  Uniform 
heat  flux,  uniform  Initial  temperature  distribution,  and  fully  developed  veloc- 
ity distribution. 
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Figure  8.  - Variation  of  x,  defined  by  t = x(tg  - t^)  + t^ , vith  x/D  for  use  with  figure 
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Figure  10.  - Concluded.  Variation  of  dimensionless  wall  and  bulk  temperatures  with  X/D  and 
Reynolds  number  for  gas  flowing  in  a tube.  Pr,  0.73.  Uniform  heat  flux,  uniform  initial 
temperature  distribution,  and  fully  developed  velocity  distribution. 


* XI.  - Variation  o'  average  Huaselt  number  defined  by  equation  (?i  ) with  X/D  and  Reynolds  number  **or  gas  flowing  in  a tube 
0.73.  Uniform  heat  flux,  uniform  initial  temperature  distribution,  and  fully  developed  velocity  distribution. 


Figure  11.  - Concluded.  Variation  of  average  NuBEelt  number  defined  by  equation  (25)  with  X/D 
and  Reynolds  number  for  gas  flowing  in  a tube.  Pr,  0.73.  Uniform  heat  flux,  uniform  initial 
tes^erature  distribution,  and  fully  developed  velocity  distribution. 
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(c)  qj'  » -0.0025  (heat  extraction  from  gas). 

Figure  14.  - Concluded.  Variation  of  Nusselt  number  with  x/D  and  Reynolds  number 
for  gas  flowing  between  parallel  flat  plates.  Pr,  0.73.  Uniform  heat  flux,  uni- 
form initial  temperature  distribution,  and  fully  developed  velocity  distribution. 


Figure  15.  - Concluded.  Variation  of  dimensionless  wall  and  bulk  temperatures  with 
X/D  and  Reynolds  number  for  gas  flowing  between  parallel  flat  plates,  Pr,  0.73. 
Uniform  initial  temperature  distribution  and  fully  developed  velocity  distribution 


Figure  16.  - Continued.  Variation  of  average  Nusselt  ruirber  defined  by  equation  (rr  ) with  X/D  ani 
Reynolds  number  for  gas  flowing  between  parallel  l'ln’  jlates.  Pr , 0.73.  Unil'orm  heat  flux,  uni- 
form initial  temperature  distribution,  and  fully  developed  velocity  distributions. 
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o r 4 6 f 10  12 

x/d 


(c)  q j ' = -O.OOPS  (heat  extract!  >ri  from  gas). 

Figure  IB.  - Concluded.  Variation  of  average  Nusselt  number  defined  by  equation 
(?.')  with  X/D  and  Reynolds  number  **or  gas  flowing  between  parallel  flat  plates. 
Pr,  0.73.  Uniform  heat  flux,  uniform  initial  temperature  distribution,  and 
fully  developed  velocity  distribution. 
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1st  loti  of  Kuaaelt  miabtr  with  X/D  and  Reynold*  nuaber  for  gas  flowing  In  a 
0.73.  ttalfcra  temperature , uniform  Initial  temperature  distribution,  fully 

- - * ■»  J J -A M V. .A.  1 m _ _ _ . . . « * 


rrcpcrt : 
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Figure  25.  - Predicted  generalized  temperature  distribution  f^r  *'l  tv  ;>f  liquid  metal  with  heat 
transfer.  Data  ■f’rom  reference  Pr,  0.01. 
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Figure  29.  - Variation  of  friction  ractor  based  'X)  static-pressure  gradient  with  X/D  and  Reynolds  number  for  flow  In  a tube 
Constant  properties. 
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Figure  32.  - Variation  of  friction  factor  based  on  static -pressure  gradient  with 
X/D  and  Reynolds  number  for  flow  between  parallel  flat  plates.  Constant 
properties. 
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